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We calculate analytically the effective mass and the quasiparticle renormalization factor in an electron liquid 
with long-range Coulomb interactions between electrons in two and three dimensions in the leading order den- 
sity expansion. We concentrate on the temperature dependence of the effective mass in the limit T/T-p «: r s •« 1 
and show that the leading temperature correction is linear in two dimensions and proportional to T 2 ln(l/T) in 
three dimensions (positive in both cases). We explicitly calculate the coefficients, which are shown to be uni- 
versal density independent parameters of the order of unity (in the high-density limit). The singular temperature 
corrections are due to the singularity in the dynamic dielectric function at &> ~ v?q and q <K 2pp. In two dimen- 
sions, we predict a non-monotonic effective mass temperature dependence and find that the maximum occurs 
at a temperature T* ~ 7>r s In -1 (l/r s ). We also study the quasiparticle renormalization factor in both three and 
two dimensions. 
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I. INTRODUCTION 

The basic postulate of Fermi liquid theory is the existence 
of a one-to-one correspondence between the states in a free 
Fermi gas and an interacting quantum system. This allows one 
to use the non-interacting language in describing quantum liq- 
uids. In particular, one can regard the interacting Fermi sys- 
tem as a gas of elementary quasiparticles. In this approach, 
the number of parameters describing the state of the system 
is less than within the exact description. Thus, an elemen- 
tary excitation is not a stationary state but a wave packet of 
stationary states which spreads with time. This leads to a 
finite life-time of elementary excitations t p away from the 
Fermi surface. However, if the inverse life-time is smaller 
than the excitation energy r p £ p » 1, one can regard the ex- 
citations as stable particles. The effective mass m* of these 
particles is renormalized by the electron-electron interactions 
and can be quite different from the non-interacting bare elec- 
tron mass (the band mass, m). The concept of the electron 
effective mass has been a subject of investigation for over 
fifty years. Surprisingly, the question of the effective mass 
temperature dependence had never been addressed until very 
recently. 12 - 3 This can be partially explained by the fact that 
most of the work was performed back in the fifties and six- 
ties, when only the three-dimensional case was of interest. In 
typical three-dimensional systems, the Fermi energy is very 
high compared to the temperatures relevant to experiments 
(i.e., in simple metals: T-p ~ 10 4 K) and thus any tempera- 
ture corrections are negligible. The Fermi energy in realistic 
semiconductor-based two-dimensional systems {e.g. Si MOS 
structures, GaAs heterostructures and quantum wells) may be 
as low as IK, which makes the issue of the temperature de- 
pendence of Fermi liquid parameters extremely important. In 
this paper, we obtain analytical results for the quasiparticle 
effective mass and the quasiparticle renormalization factor for 
two-dimensional and three-dimensional electron systems in- 
teracting via the realistic long-range Coulomb potential. Our 
results employ the standard perturbation theory expansion in 
the dynamically screened interaction and are exact in the high- 
density limit. We restrict ourselves entirely to the case of an 



ideal clean (i.e., no impurity disorder) and homogeneous elec- 
tron system with a parabolic non-interacting energy disper- 
sion. 

Before describing the main results and the structure of our 
paper, let us briefly discuss previous studies in the subject. 
The first work explicitly calculating the effective mass is due 
to Gell-Mann, 4 who derived a zero temperature correction to 
the effective mass due to the Coulomb interaction in the high- 
density limit in three dimensions. Galitski: 5 has developed a 
general scheme of calculating perturbative corrections to the 
one-particle spectrum of an interacting Fermi-system. Within 
this scheme, he derived corrections to the quasiparticle effec- 
tive mass and lifetime for the cases of both a short-ranged 
interaction and the long-ranged Coulomb interaction. These 
studies were again constrained to zero temperature. Chaplik 6 
and later Giulianni and Quinr. 7 have addressed the issue of the 
quasiparticle lifetime temperature dependence having found 
in two dimensions a non-analytic contribution to this quan- 
tity EpTp 1 cc max J^p, T 2 J In |max J^ p , rJJ. This result assures 
that the quasiparticles are well-defined excitations as long as 
T <*c Tp. Very recently, Chubukov and Maslov 1 - 2 revisited the 
problem of non-analytic corrections to the Fermi-liquid the- 
ory for the case of a short-ranged interaction. In particular, 
they showed that the leading temperature correction to the ef- 
fective mass is linear, similar to the results for the Coulomb 
interaction case, which had been reported in a short numerical 
paper earlier^ 

In this paper, we present detailed purely analytic calcula- 
tions of the effective mass renormalization by the Coulomb 
interaction. We work within first order perturbation theory in 
the screened interaction, i.e. within the random-phase approx- 
imation (RPA). We analytically derive the leading tempera- 
ture corrections to the effective mass in the low temperature 
T/Tp <k r s and high density r s <K 1 limits. In two dimen- 
sions, the leading correction is positive and linear in tempera- 
ture with the subleading term being of the order of T 2 In T and 
negative. The linear term coefficient is found to be a density 
independent universal number. In two dimensions, we pre- 
dict a non-monotonic effective mass temperature dependence. 
The point of maximum of the curve m*(T) is calculated ex- 
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plicitly and is shown to drift toward higher temperatures as r s 
increases. 

For the sake of completeness, we also calculate m*(T) 
for the case of a three-dimensional electron liquid. At high 
densities r s <k 1, the leading contribution is of the order 
of T 2 ln(l/r) and positive. As r s increases, the correction 
changes its sign and m*(T) monotonically decreases from its 
zero temperature value. 

In addition to the quasiparticle effective mass m*, another 
important many body Fermi-liquid parameter is the quasipar- 
ticle renormalization factor (the Z-factor), which is a measure 
of the quasiparticle spectral weight. In particular, the Z-factor 
defines the size of the effective Fermi surface discontinuity in 
an interacting system, and is precisely the size of the discon- 
tinuity in the momentum distribution function n{p). For the 
non-interacting Fermi gas, n(p) = (pp - p), and the discon- 
tinuity is precisely unity whereas for an interacting system this 
discontinuity is Z < 1 . Note that Z + implies the validity of 
Fermi liquid theory. To the best of our knowledge, there has 
been no consistent microscopic analytic derivation of the in- 
teraction corrections to the quasiparticle Z-factor even in three 
dimensions and zero temperature. To fill this gap we calculate 
analytically the Z-factor. The technical part of this calculation 
is found to be more complicated than the effective mass cal- 
culation. To get the correct result, one is required to use the 
exact forms of the polarizability to ensure the convergence of 
the final result. 

Our paper is structured as follows: In Sec. HI] we give a gen- 
eral introduction and derive the basic formulae for the analyt- 
ically continued self energy in first order perturbation theory 
in the screened interaction in three and two dimensions. In 
Sec. Ill CI we briefly discuss the structures of the interaction 
propagator and the polarization operator in two and three di- 
mensions. In Sec. [HI] we study the temperature dependence 
of the effective mass and also derive the asymptotic formula 
for the quasiparticle Z-factor in a three dimensional electron 
liquid. In Sec. lIVI we study the two-dimensional case and find 
that the leading temperature correction to the effective mass is 
linear and positive with the subleading term being of the order 
of T 2 In T and negative. Hence, the effective mass temperature 
dependence is non-mono tonic. We explicitly derive the tem- 
perature T* at which the maximum of m*(T) occurs. We show 
that the non-analytic contribution to the effective mass is due 
to the singularity of the polarization operator at co ~ upq and 
ij« 1. We also discuss the asymptotic behavior of the quasi- 
particle renormalization factor in two dimensions and show 
that within the RPA approximation the correction is indeed fi- 
nite and negative. We emphasize that this result is due to a 
subtle cancellation of the logarithmic singularities in the (§§)- 
derivative. These kinds of singularities persist in higher orders 
of perturbation theory. 



n. GENERAL FORMULAE 

In this section, we give basic formulae which will be used 
for actual calculations. We consider a spin degeneracy fac- 
tor of 2 throughout the paper. We denote the quasiparticle 



effective mass (bare mass) as m* (m). We use the following 
small parameters: a (3D) = e 2 /(7inv F ) and a (2D) = e 2 /(frv F ) 
in three and two dimensions respectively. These true pa- 
rameters of the asymptotic expansion are connected with the 
usual r s -parameter as follows: r^ D) = n(9n/4) 1 ^ 3 a (iD ^ and 
V2a (2D) . In what follows we will use units % = ka = 1 . 



„(2D) 



A. Renormalized spectrum 

The exact Green function for a system of interacting 
fermions can be expressed in terms of the self-energy S(e, p) 
as follows 



G (e, p) = s - £ (P) + H - S(e, p), 



(1) 



where £b(p) = p 2 /2m is the spectrum of non-interacting 
fermions and p = Ep is the bare chemical potential. The 
Green function can be rewritten as 

G(e, p) = ^ . , -, (2) 

s - £(p) + iy(s, p) 

where £(p) is the renormalized spectrum of excitations, 
y(e, p) is the quasi-particle decay rate, and Z is the residue of 
the Green function, which determines the jump in the Fermi 
distribution at p = pp. In what follows, we mostly will be 
interested in the renormalized single-particle spectrum, which 
is connected with the real part of the self-energy as follows: 



£0?) = £oCp) + Re £(£,£,) 



(3) 



where we write the self-energy as a function of a small exci- 
tation energy £ p : 

AO) P 2 PF , \ 

Sp = ~ p~ —(p-pF), 

Zm m 



£ p = E(p) - p* 



Pf 



(p - Pf) ■ 



The shift of the chemical potential of quasiparticles is deter- 
mined by the following equation 



p -p = E(0,0;//). 
From the above equations, it follows: 

^p=4 0) +Re[E( e = ^ p )-X(0,0)]. 



(4) 



(5) 



Solving the linear equation, we obtain the following formula 
for the renormalized effective mass: 

-l 



m*(T) 



d 

1 + — Re £(£,£,) 

Cttfp 



6,fp=0 



with 



1_ — ReZ(e,f p ) 

OS 



e,fp=0 



In the perturbative regime, the effective mass reads: 
m*(T) 





d d 






1 - 


ds <9£ p 


ReS(e,^) 





(6) 



(7) 



(8) 



B. Self-energy in the RPA approximation. 

In first order perturbation theory in interaction, the Matsub- 
ara self-energy can be written as (see Fig. la)£ 

S(e„, p) = -T ~ (o m ,P ~ q)£>(to m , q), (9) 

<*>m 

where s„ = n(2n + l)T is the fermion Matsubara frequency, 
co m = 2nmT is the boson Matsubara frequency, and T is the 
temperature. The function D(q, co m ) denotes the coupling to 
a collective mode (phonon, plasmon, electron-hole excitation, 
etc.), i.e., an effective interaction. 

For analytical calculations, it is more convenient to use the 
self-energy as a function of the real frequency s rather than 
the Matsubara frequency s„. It is known 9 that in some cases 
(e.g., when calculating the ground state energy in the RPA) it 
is more convenient to do the Matsubara summations first to 
avoid divergences arising form the plasmon pole. However, 
when calculating the effective mass or the quasiparticle renor- 
malization factor one can do the analytical continuation first 
and obtain finite results. In the calculation of the effective 
mass in the high-density limit, the plasmon singularity does 
not show up in the calculations at allm The calculation of the 
Z-factor is more complicated but the correct result can be ob- 
tained by keeping the exact ^-dependence in the polarizability 
(see Sec. IIII B\ . Using the standard procedure of the analytical 
continuation, one can obtain the following expression for the 
analytically continued self-energy function: 8 

d d q 



(a) 



x< Im G R (co + b, p - q) D R (-co, q) tanh 



co + e 



2T 



+ G R (co + s, p - q) Im D R (co, q) coth 



IT 



(10) 



where functions labeled with index R are retarded functions, 
i.e., functions analytical in the upper half-planes of the com- 
plex frequency and d stands for the dimensionality of space. 

Within leading order perturbation theory, one can use the 
bare electron Green function in Eq. ( 1 1 Oft , which can be written 
as 

Gf (s,p) = [s-^ + i0]-\ 

If the effective interaction is isotropic (which indeed is the 
case in the jellium model we wish to study), the integral(s) 
over the directions of q can be evaluated and we obtain the 
following expressions for the real part of the retarded self- 
energy function: 

Re E R (s, £ p ) = (s, f p ) + 2 2 (s, , 

where in three dimensions: 

+oo qi(0.) 

Zj 3D) [s, gpj = ^— j — dco Jdqq Re D R (co, q) tanh 

-oo qi (£1) 



(ID 



2T 
(12) 



+ 



+ i 



(b) 




(c) 




FIG. 1 : (a) The self-energy diagram relevant in the high-density limit 
with the thick wiggly line being the dynamically screened Coulomb 
interaction and the thin wiggly line the bare Coulomb interaction. 
The solid lines correspond to the electron Green functions; (b) The 
dynamically screened Coulomb interaction infinite series in the bare 
interaction through the bubble diagrams. A bubble corresponds to the 
noninteracting ("Lindhard") polarizability. (c) Examples of higher 
order diagrams which are negligible in the high density limit but are 
important at lower densities. 



and 



fee oo 



r(3D) 



O - q 2 /2m + upq 



X Im D R (co, q) coth 
and in two dimensions 



2T 



+oo q 2 (£1) 



Q - q 2 /2m - v^q 

(13) 

dq 



x Re D R (co, q) tanh 



co + s 



2T 



(14) 



and 



r (2D) 



dq 



V "Fl ) 



: Im D R (co, q) coth 



?,(£!) 
CO 



2T 



(15) 



In Eqs. (I12> . Jl 3i . 14l >. and dl5> . we introduced the following 
notations for the sake of brevity: 

£1 = co + s — ^ p 
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and 02 (£2) > qi (£2) > which are the solutions of the equa- 
tion Q - q 2 /2m\ = v F q. 



C. Effective interaction and the polarization operator. 

The appropriate propagator in the case of an electron liquid 
with long-range Coulombic forces between electrons is given 
by the sum of the ladder of bubble diagrams and has the typi- 
cal RPA form (see Fig. lb): 



D(o), q) 



V(q) 



1 +V(q)IL(a>, q) 



(16) 



where V(q) is the bare Coulomb interaction [V(q) = 4ne 2 /q 2 
in three dimensions and V(q) = 2ne 2 /q in two dimensions] 
and q) is the polarizability bubble: 

/dEd d x> 
— — i-G (0) ( e ,p)G (0) ( £ + W ,p + q). (17) 
(2nf ' 

At zero temperature, the polarizability was calculated by 
Lindhard 10 and Stern 11 in three and two dimensions respec- 
tively. We will need the exact expressions, which can be 
conveniently written in terms of the dimensionless parameters 
u — cjj (upq) and x = q/(2pp). In three dimensions it reads 



V (iU > fir n 

Re n (3D) («, x) = — | 1 + — [l - (jc + u) 2 \ In 
+ — [l -(x- M ) 2 lln 

Ax L J 



1 + X + U 



1 — x — u 
1 + x — u 



1 — X + u 



(18) 



and the imaginary part for u > [Imn(w, x) 
- Imfl (-u, x)]: 



Imn (3Z3) (u,x) 



it (\ \ — x\ — u) 



+ -(x-u) 2 \6(\\ +x\-u)6(u-\\ - jc|)|, (19) 

where v (3D) = mp/n 2 is the density of states at the Fermi sur- 
face. 

In two dimensions the polarizability has the following ex- 
plicit form: 



Re Tl (2D> (u, x) 



J2D) 



Re-M 



sgn (x + u) 

2x 

sgn (x - u) 



J(x + u 



2x 



■\j(x — u 



f - 1 (20) 



and the imaginary part for any u: 



Let us emphasize that both three dimensional and two 
dimensional polarizabilities are non-analytic functions at 
\u ± x\ - 1. Usually this singularity is associated with Friedel 
oscillations and Kohn-Luttinger effect, 12 i.e. with the fa- 
mous Kohn anomaly at q = 2pp and a> — > 0. In the case 
of a dense Coulomb liquid, the typical momenta are small 
x = q/(2pp) <k 1, but u = a>/(upq) can be of the order of 
unity or even larger. As we shall see, exactly this domain of 
parameters (i.e., small momentum transfers) is responsible for 
non-analytic contributions to the effective mass temperature 
dependence and the quasiparticle Z-factor. The usual Kohn 
singularity (x - 1) in the static polarizability also gives rise to 
a non-analytic temperature dependence but this effect is para- 
metrically smaller than the dynamic screening effects in the 
limit r s <K 1 . 

The issue of the temperature dependence of the polariza- 
tion operator was recently reconsidered in great details by 
Chubukov and Maslo^ (see also Refs. fl3ll and ffl ) who 
found that in the vicinity of the Kohn singularity the polariz- 
ability has a linear T correction, which is important in the case 
of a short range interaction case. In the case of the long-range 
Coulombic forces, the ^-dependence of the propagator 
becomes crucial and in the high-density limit the results are 
determined by the region x = q/(2pp) ~ r s <k 1 in two dimen- 
sions (x ~ -\ff\ in three dimensions). Hence, in the leading 
order in r s , only the region i < 1 is important. In this re- 
gion, the leading temperature correction to the polarizability 
is of the order of T 2 in any dimensionality. As we shall see 
below, the leading order temperature corrections to the effec- 
tive mass and the Z-factor are parametrically larger than T 2 . 
Therefore, in the high-density limit, the temperature correc- 
tions to the polarization bubble give negligible contributions 
to the quasiparticle spectrum and in the case T/Ep <K r s <k 1, 
one can use the zero temperature results for the polarizability 
[see Eqs. dl 8b . dl9> . ( 1201 . and ( 1211 1 in the calculations of the 
effective mass and the Z-factor. 



D. Collective modes 

The usual practice is to expand the polarizability functions 
in x «: 1 ; in which case the polarizability becomes the func- 
tion of just one variable u - a>/(vpq). The limit u <sc 1 cor- 
responds to the electron-hole branch of excitations. The cor- 
responding retarded electron-hole propagator is (in three di- 
mensions) 



+ x 3 1 



Vpq) 



(22) 



where -2pf Va (2D) is the inverse screening length. 
In two dimensions the electron-hole propagator reads 



Imn (2D) ( M ,x) JLl Rejy/l - (x + uf ~y/l - (jc - «) 2 J,(21) Df h D \co, q) = 2ne 2 q + x 2 - ix 2 



(23) 



with v (2D) - m/n being the two-dimensional density of states 
at the Fermi line. 



where x 2 — 2ppa ( - 2D) is the inverse screening length in two 
dimensions. 
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The opposite limit u » 1 corresponds to the plasmon 
branch. The spectrum of plasma waves is determined by the 
equation 

1 + V(q)Yl(a), q) = 0. 

In three dimensions the spectrum has the following well- 
known form (J^ D \q) = -^w 2 + \v^q 2 , with ojq = 4we 2 n/m. 
In two dimensions, the plasmons are gapless cJ^iq) = 



^aq + bq 2 with a - 2e 2 E F and b — 3£p/(2m). 

Within the first approximation in the interaction, the imagi- 
nary part of the polarization operator at a> » v F q is zero. This 
is exactly the source of the well-known plasmon singularity in 
self-energy calculations. However, higher order diagrams de- 
liver non-zero contributions to the imaginary part. Taking into 
account this fact, we can write down the following expression 
for the retarded plasmon propagator: 



Dpi(a),q) = 



1 



u> - co p i{q) + i0 



ll> + copi(q) 



^— — ]6(q m -q), (24) 
+ 1O 1 



where q m is the wave-vector at which the strong Landau 
damping commences. 

We shall see that the leading temperature correction to the 
effective mass both in two and three dimensions comes mostly 
from the region of u ~ 1, which is neither plasmon nor 
electron-hole region. In actual calculations, we do not sep- 
arate the screened Coulomb propagator into the electron-hole 
and plasmon branches. Moreover, for the calculation of the 
renormalization factor one is required to keep the exact polar- 
izability function (without expanding oni^ 0). 



III. THREE-DIMENSIONAL CASE 

In this section we present analytic calculations of the effec- 
tive mass and the quasiparticle Z-factor in a three-dimensional 
dense electron liquid. Throughout this section, v = mp/7T 2 
is the three-dimensional density of states and a - e 2 j {nhv-p) 
is the appropriate expansion parameter. The main results are 
Eqs. J29I and d34l > given below. 



A. Effective mass 



following expression: 
d 



ReZ(e = ^) 



1 Cdco 
J 



1 



f=o 2r cosh 2 [#] 

q 2 [oJ-E (q)] 

X J dqq Re D{aj,q) .(25) 

qi[bj-E (q)] 



It is convenient to separate the static and dynamic propaga- 
tors: 

D(oj, q) = D Bt (q) + D d Ua>, q) = D(0, q) + [D(oj, q) - D(0, q)] . 



The static propagator has the form: 

a 1 
DJx) = - ■ 



(26) 



(27) 



v x 2 - + a 

At low temperatures, frequencies a> in the integral ( I25> are 
of the order of temperature or even lower. We consider the 
following asymptotic regime of ultra-low temperatures: 

5 = to/(4E F ) ~ T/E F 

In this limit, the real part of the dynamic propagator has fol- 
lowing form [see Eq. dl 81 i. 1191 1: 



ReZ) dyn (<5,x) =- _) (x 2 + a) 



1 (x 2 + a)" + ( 



\ 2 



(28) 



Using Eq. ( 125 \ and propagators J27l > and J28l >. one can calcu- 
late corrections to the effective mass (expanding on the small 
parameter S ~ T/Ep) and see that the static contribution, in- 
deed, solely determines the renormalization of the effective 
mass at zero temperature. However, it gives temperature cor- 
rections of the order of T 2 only. The dynamic part gives zero 
contribution to the zero temperature effective mass, but gives 
parametrically larger temperature dependent correction. The 
final result reads: 



m* (T) - m a 1 n ( T \ Ev 
= — ln- + — — In — . 

m 2 a 96\E F T 



(29) 



The first term indeed coincides with the old result of Gell- 
Mann. The non-analytic temperature correction given in the 
second term of Eq. d29i is a new result. Let us emphasize 
that the leading correction is positive only in the high-density 
limit. In Ref. [3], it was shown that the leading T 2 In T correc- 
tion changes sign at lower densities (a (3£)) ~ 1). The density 
dependence comes from large momentum transfers q, in par- 
ticular from the vicinity of the 2/?F-anomaly, which becomes 
increasingly important at lower densities. 



In the high density limit, the correction to the effective mass 
is determined by the on-shell equation which means that 
we can put s = g p and study the self-energy as a function of 
just one variable [see Eqs. 1181 1. On the shell, we have the 



B. Z-factor 

Unlike in the effective mass calculation of the preceding 
section, the quasiparticle Z-factor can not be calculated within 
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the on-shell method, and therefore the problem is more com- 
plicated. In particular, one has to consider both contributions 
to the self-energy given in ill 2b and Jl 31 . The first contribution 
at zero temperature reads: 



T-Ei(e,f = 0) 
os 



OO CO 

= — I dxx t duReD (u, 



-v) 



x— [0(l-|«-*l)sgnii]. (30) 

In the three-dimensional case, the integral over u reduces to a 
(5-function integration and we have: 



T-Zi(c,f = 0) 
os 



-I 

-o v: 



D(0,jc)- ^ ReD(l +x,jc) 



--ReZ)(l-x,Jc) 



(31) 



From Eq. d 1 3I >. we derive the second contribution: 



— I2(e,f = 0) 



OO OO 

= — I dxx I lmD(u,x) 



1 - (u - xf 1 - (u + x) 2 



(32) 



Using Eqs. d 1 8t > . ( I19> , and (I24l i. one can evaluate integrals (13 1 i 
and d32t and obtain 



= 0) 

os 



E=0 



OO 

I J 1 + M 2 



In 



1 



1 - m atan i 



(33) 



Evaluating the remaining integral numerically, we derive the 
final result for the quasiparticle renormalization factor: 



Z = 1 - 1.067a. 



(34) 



This asymptotic result is in a very good agreement with nu- 
merical simulations. 9 



Sec. HII B> . The only difference is the square root function, 
which appears in all two-dimensional expressions [see, e.g., 
Eqs. (I14> and dl5H . This square root singularity is identical 
to the singularity in the Stern's polarizability function given 
in Eqs. J20b and i21i . The combination of these two singular- 
ities leads to a stronger temperature dependence (linear as we 
shall see) as compared to the three-dimensional case. From 
Eqs. © and dl4> . we get the following on-shell expression 
for the two-dimensional electron effective mass: 



m* - m v £ doj 1 
m 2n J 2T C osh 2 



-/(«), 



(35) 



IT 



where we have introduced the following integral: 



I 

xi(eo) 



dx- 



x Re D (cj, x) 



^-x\(u)\ [x\{a>)-xi\ 



(36) 



where x — qj (2p F ) and are the solutions of the equation 

- xj = 1. In the limit of low temperatures: X\ ~ \6\ = 

to\/(4E F ) and;c 2 ~ 1. 

Again we re-write the propagator as a sum of static and 
dynamic terms ( I26l l. The two-dimensional static propagator 
has the form 



D st (x) = — 



v x + a 



(37) 



The real part of the dynamic propagator in the limit | a>\ ~ T <k 
E F reads: 



ReD dyn (6,x) =-- (almn) 2 (x + a)~ l 
v 



x 



[(j + a) 2 + (almn) 2 ] ' , 



(38) 
(39) 



where Im II is defined by Eq. Mil . 

Expanding in the small parameter 5, we obtain the contri- 
bution to the integral \36\ due to the static propagator 



a 

us) = - 

v 



in I 

a 



— t In — 

2a 2 \6\ 



and the "dynamic part'' 



IV. TWO-DIMENSIONAL CASE 

In this section we present analytic calculations of the ef- 
fective mass and quasiparticle Z-factor in a two-dimensional 
dense electron liquid. Throughout this section, v = mjn is 
the two dimensional density of states and a = e 2 / (Hof) <k 1 
is the appropriate expansion parameter. The main results are 
Eqs. d40l . lEffl. and d46l given below. 



A. Effective mass 

The calculation of the effective mass in two dimensions 
is analogous to the calculation in three dimensions (see 



hyn(S) = - 
V 



_n\6\_ j_ 



2 a 



\5\ 



Evaluating the elementary integral, we obtain the final result 
for the temperature dependent effective mass in the second 
leading order in temperature: 

^ = -! h i^(il^(ir^, ( 40) 

m 7i a 4 \E F ) 48a\E F ) T 

Let us emphasize that Eq. (I40l > is valid in the low temperature 
and high-density limit: T/Ep <K a <K 1 and only for sub- 
thermal particles: s <k T. We see that the leading term is 
linear in temperature and the coefficient is a universal density 
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independent number. This universal behavior is true only in 
the high density limit. There are other linear- T contributions, 
such as the one due to the temperature dependence of the po- 
larizability in the vicinity of the Kohn singularity (considered 
in the paper of Chubukov and Maslov 2 for short-range inter- 
actions). In the case of the long-range Coulomb interaction, 
the Kohn anomaly leads to a linear- T term proportional to the 
Coulomb expansion parameter a. Similar a-dependence of 
the linear slope was discovered in RPA numerical calculations 
in Ref. In the high-density limit, this density dependent 
linear- T term is asymptotically smaller than the main univer- 
sal contribution [the second term in Eq. J40H and therefore 
not shown in Eq. J40I . 



m*(T)-m*(0) 
m 




FIG. 2: Asymptotic effective mass temperature dependence [see 
Eq. <40H for three different densities n\ > n 2 > n 3 . The slope of 
the curves is universal at high densities [ -\/n/(me 2 ) » 1]. 

From Eq. (I40t . we see that the effective mass temperature 
dependence is non-mono tonic. A maximum occurs at a tem- 
perature T*, which within the logarithmic accuracy has the 
form: 



T" 



6ln2 a 
5tt In 



0.26 



ln± 



<sc a. 



This result is formally within the limits of applicability of our 
theory. We see that the point of maximum of the curve m*(T) 
drifts toward higher temperatures as the density decreases. 
This tendency is preserved at lower densities as well (within 
the RPA approach). Such a maximum in m*(T) and a density 
dependent T* were also discovered in our recent numerical 
calculation. 3 

B. Z-factor 

The analytical calculation of the Z-factor in two dimensions 
is technically a very demanding problem. The mixture of two 
singularities, the Kohn singularity in the polarizability and the 
identical square root singularity in Eqs. 14l > and dl5l > arising 
from the two-dimensional phase space, leads to a complicated 
structure of the integrals in Eqs. H41 and d!5l >. each being 
a truly divergent quantity. The logarithmic divergence gets 
cancelled (at least within the RPA), but to see this cancella- 
tion one is required to keep the exact x and u dependences 



in the Stern's polarizability function. Moreover, the technical 
method used in the three-dimensional calculation of the Z- 
factor [see Sec. II 11 Bl Eq. O0H is not applicable here because 
of the square-root singularity. 

Let us now study the two-dimensional Z-factor in more de- 
tails. The quasiparticle renormalization factor is determined 
by the energy derivative of the self-energy. The latter can be 
written as a sum of two terms [see Eqs. fl!4i and dl5H : 



(Cf = 0) 

os 



£^0 



oo oo 

I dxx I 

2n J J 



du 



— Re D (u, x) 
du 



HI - \u-x\) + 6(1 - \u + x\) 



yi - (u - xj 1 ^\ - {u + xj 1 



(42) 



and 



ds 



£ 2 (e,£ = 0) 



oo oo 



G(\u-x\ - 1) 



du 



du 



Im D (u, x) 



]u + x\ - 1) 



V(« - x) 2 - 1 V(« + x) 2 - 1 



■(43) 



The frequency dependence (hence, the M-dependence) of the 
propagator D(u, x) is due to the polarizability II(m, x), which 
contains exactly the same square root functions as the ones 
in Eqs. J42> and J43i . This leads to a logarithmic divergence 
of each of the above integrals at u — 1 ± x. The "singular" 
contributions have the following forms: 



os 



£^0 



— J" dxx J" 



du 



o o 
dlmf! 



[Re D (u, x)f 
Vl - (u - x) 2 



du 



■Imfl 



x+ - Ren 

v 



(44) 



(41) and 



5 fi 



(sing) 



= 0) 



e^0 



oo u 

— J" dxx J" 

l-.t 

5Ren. 



[Re D (u, x)] 



V(« 



du 



Imfl 



a 



1 



x+ - Ren 

v 



■ (45) 



Each of these quantities is logarithmically divergent 
(-^r) ~ +alns — > co. We emphasize that in two dimensions 
the real and imaginary parts of the polarizability have almost 
identical analytic structures, in contrast to the three dimen- 
sional case in which the imaginary and real parts are basically 
independent functions with quite different properties. Using 
Eqs. 1441 and J45> . one can check that this "symmetry" of the 
two-dimensional polarizability leads to an exact cancellation 
of the logarithmic divergence and to a finite result. 

Let us emphasize that this kind of dangerous singularities 
appear in any order of the perturbation theory in interaction 
(see, e.g., Fig. lc). It is not a priori obvious how (and if) 
the singularity, which is cut-off only by temperature or energy 
e, is cancelled in higher order diagrams. We do not have a 
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general argument for why the divergence must cancel in each 
order, but we do know that they cancel to this order. It is 
essential to clarify this point to assure that the quasiparticle 
Z-factor does not vanish logarithmically Z~'(e) ~ lne and to 
make certain that the usual Landau Fermi liquid theory is pre- 
served in two dimensions. This issue is currently being stud- 
ied by us. We believe that the Fermi liquid theory is preserved 
but it needs to be demonstrated explicitly. 

Within the RPA, the zero temperature Z-factor can be 
proven to be finite. 15 One can formally define the quasiparticle 
Z-factor at finite temperatures via relation Q. Studying the 
leading temperature correction to the energy derivative of the 
self-energy is quite similar to the case of the on-shell deriva- 
tive case. The leading term can be shown to be linear and 
negatived 



Z(a, T) sal — 




where c is a constant of the order of unity. 

V. CONCLUSION 

In this work we have developed the analytic leading order 
theory for the temperature dependent quasiparticle effective 
mass, m*(T), and the quasiparticle renormalization factor for 
two- and three-dimensional interacting electron systems. Our 
results are asymptotically exact in the low temperature high- 
density limits for the case of the realistic long-range Coulomb 
interaction, and thus we are complementary to the recent an- 
alytical work of Ref. 01, which considers a short-range re- 
pulsive interaction. It is interesting to note that m*(T) has an 
unexpected linear- T correction (rather than T 2 ) both in our 
theory and in the theory of Chubukov and Maslov; but in our 
case the correction is positive opposite to the short range case. 
This immediately leads to the conclusion that the leading cor- 
rection to Cy/T (where Cy is the specific heat) is not univer- 
sal — our long range interaction produces a positive linear- T 
term in the leading order in contrast to the negative sign ob- 
tained in Ref. |2|]. This unexpected linear- T term appears due 
to the non-analyticity of the polarizability function. This non- 



analyticity has potentially important consequences for quan- 
tum critical phenomena as discussed recently in Ref. 1 17]. 

Our analytic results are also in agreement with recent nu- 
merical studies of the temperature dependent effective mass, 3 
which employed the random-phase approximation at lower 
densities. It was shown that the linear- T correction persists 
at lower densities as well (at least within the RPA) and the 
qualitative behavior of the effective mass remains the same 
with the only difference being the density dependent slope of 
the m*(T) curve at T — > and r s > 1 (for high densities, 
it was shown to be density independent in agreement with 
the analytical results reported in the present paper). As the 
RPA is believed to be qualitatively reliable at lower densities 
as well, we expect our results to be quite general and quali- 
tatively applicable to realistic two-dimensional electron sys- 
tems. We also should point out that although the temperature 
dependent effective mass renormalization has only been cal- 
culated numerically very recently, 3 there is a vast literature of 
numerical studies of zero temperature many-body effects in 
the three-dimensional and two-dimensional interacting elec- 
tron systems. We cite in this context only two rather compre- 
hensive references: Ref. 1 18] for three-dimensional systems 
and Ref. 1 15] for two-dimensional systems. 

We have proved that the quasiparticle Z-factor in two di- 
mensions is finite at least within the random-phase approxi- 
mation. However, we would like to emphasize that each order 
in perturbation theory does contain a dangerous logarithmic 
singularity in this quantity. It is known that in two dimensions 
higher order diagrams may contain important effects (see, e.g. , 
Ref. 11911 . in which it was shown that the static Kohn-Luttinger 
effect 12 in two dimensions is "hidden" in third order perturba- 
tion theory only). It is therefore essential to prove that the 
cancellation of the logarithmic singularity (which would oth- 
erwise lead to a logarithmically vanishing Z-factor and to a 
marginal Fermi liquid 20 ) takes place in higher orders. This 
important question will be considered elsewhere. 
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